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Abstract — This paper investigates the performance of an or- 
thogonal frequency-division multiplexing (OFDM)-based cogni- 
tive radio (CR) spectrum sharing communication system that 
assumes random allocation and absence of the primary user's 
(PU) channel occupation information, i.e., no spectrum sensing 
is employed to acquire information about the availability of 
unused subcarriers. In case of a single secondary user (SU) 
in the secondary network, due to the lack of information of 
PUs' activities, the SU randomly allocates the subcarriers of 
the primary network and collide with the PUs' subcarriers 
with a certain probability. To maintain the quality of service 
(QoS) requirement of PUs, the interference that SU causes onto 
PUs is controlled by adjusting SU's transmit power below a 
predefined threshold, referred to as interference temperature. In 
this work, the average capacity of SU with subcarrier collisions 
is employed as performance measure to investigate the proposed 
random allocation scheme for both general and Rayleigh channel 
fading models. Bounds and scaling laws of average capacity with 
respect to the number of SU's, PUs' and available subcarriers 
are derived. In addition, in the presence of multiple SUs, 
the multiuser diversity gain of SUs assuming an opportunistic 
scheduling is also investigated. To avoid the interference at the 
SUs that might be caused by the random allocation scheme and 
obtain the maximum sum rate for SUs based on the available 
subcarriers, an efficient centralized sequential algorithm based on 
the opportunistic scheduling and random allocation (utilization) 
methods is proposed to ensure the orthogonality of assigned 
subcarriers. 

Index Terms — Random allocation, subcarrier collision, OFDM- 
based cognitive radio, centralized scheduling, spectrum sharing, 
capacity, multiuser diversity. 



I. Introduction 

ADVANCES in communications technologies entail de- 
mands for higher data rates. One of the popular solutions 
to fulfill this requirement was to allocate additional bandwidth, 
which unfortunately is not anymore viable due to spectrum 
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scarcity. Recent spectrum measurement campaigns, performed 
by agencies such as Federal Communications Commission 
(FCC), reported that the radio-frequency (RF) spectrum is 
being used inefficiently (see e.g., |l) and references cited 
therein). Therefore, the idea of cognitive radios (CRs) was 
advanced as a promising approach for the efficient utilization 
of RF spectrum [2|. Generally, in CR networks the usage of 
spectrum by cognitive (secondary) users is maintained by three 
approaches. In interweave cognitive networks, primary and 
secondary users are not allowed to operate simultaneously, 
i.e., the secondary user (SU) accesses the spectrum while the 
primary user (PU) is idle. In underlay cognitive (spectrum 
sharing) networks, PUs are allocated a higher priority to use 
the spectrum than SUs, and the coexistence of primary and 
secondary users is allowed under the PU's predefined interfer- 
ence constraint [3]-[10], also termed interference temperature. 
In overlay cognitive networks, SUs and PUs are allowed 
to transmit concurrently with the help of advanced coding 
techniques | [TT| . 

One of the most challenging issues in the implementation 
of CR networks is to know whether at a certain physical 
location and moment of time the RF spectrum is occupied 
by PU(s), i.e., if there is a sensing mechanism in place for 



the available spectrum [12], [13]. The challenge in deploying 
such a spectrum sensing mechanism is due to the uncertainties 
ranging from channel randomness at device and network- 
level uncertainties, to the hidden PU problem and sensing 
duration [74). fl5) . There have been numerous studies to deal 
with these issues. In fl6) and references therein, a compact 
survey of the spectrum sensing algorithms and CR applications 
along with the design and implementation challenges are 
classified properly. 

To understand the performance limits of a spectrum sharing 
system, SU capacity is a very useful performance measure. 
The ergodic and outage capacities of CR spectrum sharing 
systems in Rayleigh fading environments are studied in JS), 
and a comprehensive analysis considering various combina- 
tions of power constraints under different types of channel 
fading models is performed in (9). In fl7) , considering a 
point-to-point communication scenario, the expressions for the 
average capacity of a single SU assuming the existence of a 
single PU and no PU's interference are derived for different 
channel fading models such as Rayleigh, Nakagami-m and 
Log-normal. As an extension of |17|, in the SU capacity 
assuming PU's interference with imperfect channel knowl- 
edge, and the average bit error rate over Rayleigh channel 
fading were derived. The ergodic sum capacity of CRs (SUs) 
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with multiple access and broadcast fading channels with long- 
term average and short-term power constraints was established 
using optimal power allocation schemes in [3|. Opportunistic 
SUs scheduling yields multiuser diversity gain due to the 
channel fading randomness, and it has been well studied 
in conventional wireless systems |18|, fl9). The multiuser 



diversity analysis was conducted for spectrum sharing systems 
in (5), pO) , and interweave CR networks (see (TT) and the 
references cited therein). 

In orthogonal frequency-division multiplexing (OFDM)- 
based systems, frequency spectrum is a precious and scarce 
resource that is divided into non-overlapping bands, called 
subcarriers, and which are assigned to different cells and/or 
users. Starting with the early deployment of cellular mobile 
communication networks, efficient sharing of the available 
radio spectrum among the users has represented an important 
design problem. In conventional OFDM(A)-based systems, 
stochastic subcarrier collision models have been proposed 
to investigate the performance of various scheduling and 
deployment methods, and to assess the inter-cell-interference 
(ICI) for cell-edge users ||20)-j22). 

Considering the challenges and implementation issues in 
CR networks in terms of spectrum sensing and subcarriers 
scheduling, the existing studies motivate us to investigate the 
performance of a primitive (basic) OFDM-based CR system in 
which the SUs randomly (blindly) utilize the available subcar- 
riers assuming that some of the subcarriers are utilized by the 
PUs. This paper focuses on such a communication scenario 
that assumes random allocation and no spectrum sensing. 
An immediate challenge to be addressed is the fact that the 
SU's subcarriers collide with PUs' subcarriers. However, there 
are no studies available to assess the effect of subcarrier 
collisions in such CR spectrum sharing systems. Therefore, 
the requirement for a more comprehensive system analysis 
including the development of a stochastic model to capture the 
subcarrier collisions and protection of the operation of PUs in 
an OFDM-based CR spectrum sharing system turns out to be 
crucial. In the presence of multiple SUs, due to the random 
subcarrier allocation scheme, collisions will occur among the 
subcarriers used by the SUs in addition to the collisions 
with the subcarriers used by the PUs. The collisions among 
the SUs' subcarriers will decrease the system performance 
drastically. To overcome this issue, this paper presents also 
an efficient centralized algorithm that sequentially assigns the 
randomly selected subcarrier sets to the SUs while maintaining 
the orthogonality among these sets, to avoid collisions between 
their subcarriers. In the proposed centralized algorithm, the 
opportunistic scheduling of users, which yields multiuser 
diversity gain, is employed and the performance limits of the 
system in terms of multiuser diversity gain and sum capacity 
of multiple SUs are studied. 

The main results of this paper are next summarized: 
• A random subcarrier allocation method, where an arbi- 
trary 771th SU randomly utilizes i 7 ^ subcarriers from an 
available set of F subcarriers in the primary network, 
in an OFDM-based system is proposed. In the proposed 
scheme, the SUs do not have knowledge about the 
PUs' subcarriers utilization, i.e., no spectrum sensing is 



performed. Therefore, with some probability collisions 
between the subcarrier sets of PUs and SU occur. It 
is shown that the subcarrier collision model follows a 
multivariate hypergeometric distribution. 

• Considering the average capacity as performance mea- 
sure, the SU average capacity expressions under the 
interference constraint of PUs in the case of single or 
multiple PU(s) are derived. Upper and lower bounds on 
average capacity are derived. It is found that the average 
capacity of the rath SU scales with respect to the number 
of subcarriers in the sets F, F,f and F% as[j]9 (1 + 1/F), 
6 (l — F,f ) and O (i 7 !^), respectively. Furthermore, the 
convergence rate of average capacity as F goes to infinity 
is found to be logarithmic. 

• To find the probability density function (PDF) and outage 
probability (cumulative distribution function (CDF)) of 
the SU capacity, which is the sum capacities of subcarri- 
ers with "interference" and "no-interference" from PU(s), 
the characteristic function (CF) and moment generation 
function (MGF) approaches are in general used to obtain 
the PDF and CDF of sum of variates [23|. However, the 
obtained PDF and CDF for the capacity of the 7th sub- 
carrier for "interference" and "no-interference" cases are 
too complicated and intractable using the aforementioned 
approaches. Therefore, by using the moment matching 
method, the PDF and CDF of the 7th subcarrier capacity 
are approximated by a more tractable distribution, namely 
the Gamma distribution. There are various reasons for 
using the Gamma approximation such as being a Type- 
Ill Pearson distribution, widely used in fitting positive 
random variables (RVs), and its skewness and tail are 
determined by its mean and variance |24|-|26|. Even 
though the Gamma distribution approximation makes the 
analysis much easier to track the sum of capacities of all 
collided and collision-free subcarriers, we end up with 
a sum of Gamma variates with some of the shape and 
scale parameters equal or non-equal, and not necessarily 
integer-valued. This constraint stems from the fact that 
individual PUs can have distinct or the same transmit 
power for their subcarriers. In such a case, there are 
no closed-form expressions for the PDF and CDF of 
SU capacity. Fortunately, Moschopoulos |27| in 1985, 
proposed a single Gamma series representation for a sum 
of Gamma RVs with the scale and shape parameters 
having the properties mentioned above. Utilizing this nice 
feature of Moschopoulos PDF, the PDF and CDF of SU 
capacity are obtained. 

• Using extreme value theory, the asymptotic analysis of 
multiuser diversity is investigated. The analysis con- 
ducted at this stage reveals a novel result: the limiting 
CDF distribution of the maximum of 1Z RVs following 
a common Moschopoulos PDF and CDF converges to a 
Gumbel-type extreme value distribution as 1Z converges 
to infinity. 

1 Where F stands for the total number of available subcarriers in the primary 
network, and _F,f and are the number of subcarriers of the nth PU and 
the mth SU, respectively. The notation 0(-) is introduced in Definition^ 
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• A centralized sequential algorithm based on random 
allocation (utilization) and assuming an opportunistic 
scheduling method is proposed for scheduling the subcar- 
riers of multiple SUs while maintaining their orthogonal- 
ity. The probability mass function (PMF) of the number 
of subcarrier collisions for the mth scheduled SU in the 
algorithm is derived. In addition, the proposed algorithm 
is compared with the case, where the SUs are selected 
arbitrarily, i.e., no multiuser diversity gain is exploited. 
Last but not least, to present the impact of collisions 
among the SUs' subcarriers on the sum capacity of 
SUs, simulation results are provided and compared with 
the centralized algorithm performance with and without 
opportunistic scheduling. 

The rest of the paper is structured as follows. In Section [IT] 
some essential mathematical preliminaries and definitions are 



by 



provided. The system model is presented in Section III 



The SU capacity analysis over arbitrary and Rayleigh fading 



channels is investigated in Section IV The multiuser diversity 
gain in the opportunistic scheduling of SUs is studied in 
Section [V] Section VI presents a centralized algorithm for 
orthogonal subcarrier scheduling of SUs. The numerical and 
simulation results are given in Section |VTI] Finally, concluding 
remarks are drawn in Section fVlIII 



II. Mathematical Preliminaries and Definitions 

In this section, the hypergeometric distribution and some 
important definitions that are frequently used throughout the 
paper are provided. 

Proposition 1 (PMF of Number of Subcarrier Collisions): 
When the mth SU randomly utilizes (allocates) subcarriers 
from a set of F available subcarriers without replacement, 
and Ff subcarriers are being used by the nth PU, then the 
PMF of the number of subcarrier collisions, k nm , follows the 
hypergeometric distribution, k nm ~ HYPG(F^, Ff, F), and 
is expressed as: 



F 



Fi 



F s \k \F S - k 



F-Fl 



where the notation ( ) stands for the binomial coefficient. 
The average number of subcarrier collisions is 



m n 

F 



where E [•] denotes the expectation operator. 

Proof: The proof can be readily shown by interpreting 
the process of allocating the subcarriers as selecting balls 
from an urn without replacement. Furthermore, the expected 
value of the number of subcarriers is obtained from E [k nm ] = 
k nrn p{k nrn ). H 
In the case of multiple PUs, the mth SU might have 
subcarrier collisions with up to N PUs. Let k m = 
[k lm , k 2m , ■ ■ ■ , k^ m , kf m ] T £ ^of" 1 represent the number of 
collisions of the mth SU with TV PUs and with the collision- 
free subcarriers, kf m . Then, the (joint) PMF of k m is given 



Pr(K m = k m ) = p(k m ) = 



Ff 



Ff 



F f \(F 



Ff 

fcjVm 
-1 N 

n 



p _ y^JV f p 
F p 



F 

F s 



(1) 



f = r — 2~2n=i Fn stands for the number of 



where Ff = F 
free subcarriers in the primary network. One can observe 
that k m follows a modified multivariate hypergeometric dis- 
tribution k m — M-HYPG(i^, F p , F), where F p = 

[Ff, Ff, ...,Ff,F f ] T £ and the support of k m is 

given by: 



N 



kt — F b 



lim 



and k nm £ (F% + Ff - F) + , . . . , min {F%, Ff} 



where (x) + = max{0,x}. 

Definition 1 (Rate of Convergence R28\l): An infinite se- 
quence {An] converging to the limit A is said to be loga- 
rithmically convergent if 

|AA n+1 | \A n+1 -A\ 

. . , . and lim r— , 

|AA n | n^vo \A n -A\ 

both exist and are equal to unity, where AA n = A n +i — A n . If 
only linin^oo |AA n+1 | / |AA n | = 1 holds, then the sequence 
{A n } converges sublinearly to A. 

Definition 2 (Knuth's notations fl29| /J: Let f(n) and g(n) 
be nonnegative functions. The notation: 

• f(n) — 0(g(n)) means that there exist positive constants 
c and no such that f(n) < cg(n) for all n > uq. 

• f{n) — fl(g(n)) means that there exist positive constants 
c and no sucn that f(n) > cg(n) for all n > no, i.e., 
g(n) = 0(f(n)). 

• /(n) = 8(<?(n)) means that there exist positive constants 
c, d and no such that cg(n) < f(n) < c'g(n) for all 
n > no, i.e., both f(n) = 0(g(n)) and f(n) = Q,(g(n)) 
hold. 

Definition 3: The capacity of mth SU with Ff t subcarriers 
is defined as the summation of capacities for each subcarrier. 
Let 5 m j be the signal-to-interference plus noise ratio (SINR) 
for the ith subcarrier of the rnth user, then the SU capacity is 
given byj^J 



C m = J2\og{l + S m ^). 

i=l 



and v 



Definition 4 (Capacity with Collisions): Let S, 
be the SINR for the ith subcarrier of the mth SU with "inter- 
ference" and "no-interference" from the nth PU, respectively^] 

2 A11 logarithms in the following are with respect to the base e unless 
otherwise stated. 

^^m, 1 ! ' s m£ leed the signal-to-noise ratio (SNR) for the ith subcarrier. 
However, to emphasize the subcarrier collision and collision-free cases, it 
is called SINR with "no-interference" from PU throughout the paper. 
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Fig. 1, M SUs transmit to the secondary base station (SBS) using the subcarriers in the primary network with subcarrier collisions following the hypergeometric 
distribution for accessing PUs' subcarriers, [(- -): Interference-link (channel), (-): Desired-link (channel)]. 



If k nm subcarriers of the mth SU collide with the nth PU's 
subcarriers, then the capacity of SU in Definition [3] with 
subcarrier collisions can be redefined as 



C 



^ = glog(l + ^)+^log(l + ^), 



where k nm and kf m = F^ — k nm are hypergeometric RVs that 
denote the number of collided (i.e., interference) and collision- 
free (i.e., no-interference) subcarriers between the nth PU 
and the mth SU, respectively. The superscript "1" indicates 
that collisions occur with only single PU's subcarriers (any 
arbitrary nth PU) in the primary network. The SU capacity 
expression in case of multiple N PUs is given in ( |12) . 

III. System Model 

The system model is illustrated in Figure [T] where the 
primary and cognitive (secondary) networks consist of N 
PUs with a primary base station (PBS) and M SUs with a 
secondary base station (SBS), respectively. To preserve the 
quality of service (QoS) requirements of PUs in a spectrum 
sharing communication network, the interference power levels 
caused by the SU-transmitters at the primary receiver (PBS) 
must not be larger than a predefined value (\P \, i = 1, . . . , F) 
for each subcarrier, referred to as the interference tempera- 
ture (IT). It is assumed that there is no correlation among 
the subcarriers. Nonetheless, due to the inherent nature of 
random allocation (utilization) method and the high number 
of available subcarriers in practice, the probability of a SU to 
select consecutive subcarriers, which are practically correlated, 
would be considerably negligible. 

The channel power gains from the mth SU to SBS and 
PBS are denoted by h m and h mp , respectively. Similarly, g n 
and g ns represent the channel power gains from the nth PU to 
PBS and SBS, respectively. All the channel gains are assumed 
to be unit mean independent and identically distributed (i.i.d.) 
flat Rayleigh fading channels. The channel power gains are 
hence exponentially distributed with unit mean. Further, to 
have a tractable theoretical analysis, it is assumed that perfect 



information about the interference channel power gains, h mp , 
is available at SUs. The SUs can obtain this information, 
referred to as channel side information (CSI), through various 
ways, e.g., from the channel reciprocity conditioijj j6), fio) , 
or from an entity called mediate band or CR network manager 
between the PBS and SU Q. The thermal additive white 
Gaussian noise (AWGN) at both PUs and SUs is assumed 
to have circularly symmetric complex Gaussian distribution 
with zero mean and variance rj, i.e., CA/"(0,n). Throughout 



the paper, the parameters h r , 



g n:i and g ns . t denote 



the channel power gains associated with the ith subcarrier. 
Furthermore, for the sake of analysis simplicity, the value 
of IT is assumed to be the same for all subcarriers in the 
system and available at the SUs, and the transmit power of 
each user (either PU or SU) is the same for all its subcarriers, 

I.e., Pn.i — Pn &nd Prn,i Prrv 

The total number of available subcarriers in the primary 
network is denoted by F. The subcarrier set of each PU is 
assumed to be assigned by preserving the orthogonality among 
the sets of subcarriers for all PUs, F.^ for n = 1, . . . , N. SU 
randomly allocates the subcarriers from the available subcar- 
riers set F without having access to the information about the 
channel occupied by PUs. Therefore, SU will collide with the 
subcarriers of the PUs with a certain probability. Subcarrier 
collisions occur when SUs employ subcarriers which are in 
use by PUs, and the probabilistic model for the number of 
subcarrier collisions follows a multivariate hypergeometric 
distribution. 

it is 



During the evaluation of SU capacity in Section IV 



assumed that there is only a single SU (any arbitrary mth SU) 
in the cognitive network, and the collisions occur between the 
subcarriers of the SU and PUs due to the random allocation 
scheme. This set-up can also be easily extended to multiple 
SUs with the assumption of no mutual interference among 
SUs. However, such a framework would not be practical, 

4 With the assumptions of channel reciprocity and pre-knowledge of the 
PBS transmit power level, SU can estimate the received signal power from 
PBS when it transmits j6j. 
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since due to the random allocation method, the likelihood of 
the same subcarriers being allocated to multiple SUs will be 
quite high. To avoid such a scenario, an efficient allocation 
of SUs' subcarriers is needed to preserve the orthogonality 
among SUs subcarriers. Therefore, an centralized algorithm, 
which sequentially allocates the subcarriers to multiple SUs 
based on the random allocation method, while maintaining or- 
thogonality among SUs' subcarriers, is proposed and analyzed 
in Section [VI] 



IV. Capacity of Secondary User 

In this section, the average capacity of a single SU including 
the bounds and scaling laws with respect to the number of 
subcarriers for the case of an arbitrary channel fading model is 
investigated. Then, the Rayleigh channel fading model is used 
to study the impacts of the system parameters and to evaluate 
the expressions for the PDF and CDF of SU capacity. 



A. Analysis of SU Average Capacity for General Fading 

Theorem 1: The average capacity of the mth SU in the 
presence of a single (nth) PU is given by 



F s f 

77? 

~F 



F„ ( E 



a 1. 



E [C, 



NI- 

m.i 



FE [C, 



Nil 

m,zl 



where variables C 1 ^ and represent the ith subcarrier ca- 
pacity of the mth SU with "interference" and "no-interference" 
from the nth PU, respectively. In the case of Rayleigh channel 



a 



I.n 



and E [C^] are given in (13) and {H), 



fading, E 
respectively. 

Proof: The proof is given in Appendix [A] ■ 
Corollary 1: The average capacity of mth SU in the pres- 
ence of N PUs is given by 



EfC7„ 



F s 

m 



N 



5>„ P E C% +F f E[C^} 



Proof: The proof is given in Appendix [B] ■ 
1 ) Bounds on the Average Capacity: In this section, certain 
bounds on the average capacity of SU will be established. 
Intuitively, representing the relation of order between the 
average capacity of the ith subcarrier with PU's "interference" 
and "no-interference" as E C^™ < E [C^] , the naive 
upper and lower bounds on the SU average capacity can be 
expressed as 



FZE 



C 



I.n 



< 



E [GJJ < F*E [C™] 



(2) 



which states that the upper bound, in the best case, is when 
all SU's subcarriers are collision-free, i.e., all subcarriers are 
interference-free, kf m — F^. Similarly for the lower bound, 
all SU's subcarriers are colliding with the PU's subcarriers, 
i e k — F s 

However, the maximum and minimum number of subcarrier 
collisions might not be necessarily F^ and 0, respectively. The 
following general result holds. 



Corollary 2: Tight upper and lower bounds on the average 
capacity of SU in the presence of a single PU are given by: 



.maxm 



a 



I.n 



kf^E [C»% < E [Cl] 



< 



k mm ¥ 



c. 



kf^E [C» r 



J m,i\ 



where /c™^f and fc™'" represents the maximum and minimum 
number of subcarrier collisions, respectively, and are defined 

as k™ = {F,l + FP-F) + and fc™T = min{F^,F^}. 

A Ion WH IpS jumin „ nr \ l.min pS j,max 

aiso, ny m — r m K nm anu Ky m — r m n nm . 

Proof: The number of subcarrier collisions does 
not depend only on SU's subcarriers but also on 
PU's subcarriers. Therefore, the support region of 
k nm , considering the PU's subcarriers, is given by 
{(F^ + Fn -F) + ,..., mm {Ff n ,FP}}. Using this 
support region, the bounds are established. ■ 
It is worth to note that the naive upper bound, given in pi, 
on the average capacity is the limit point of capacity as the 
number of available subcarriers F goes to infinity. Formally, 



Km E[<&] =FgE [a 



Nil 
m : i\ 5 



which states that for a fixed number of PU's subcarriers as the 
number of available subcarriers increases, the average capacity 
converges to the case where no SU's subcarrier collides. 

2) Scaling Laws for the Average Capacity: 

Corollary 3: The average capacity of the mth SU in the 
presence of a single PU scales with respect to the number 
of subcarriers F, Fg and F£ as 9 (1 + 1/F), 9 (F%) and 
9 (l — F^), respectively. 

Proof: Using the Knuth's notation from Definition [2] one 
can infer that 



lim 

F-s-oo 1 



E[C4 



1_ 

F 



lim ■ 

F->oo 



F~ 



( 



E 



C 



I.n 



E[C, 



NI] 

rn,i\ 



■ F%E 



F^E [CZ] > 0. 



F 



Following the same approach, one can establish the scaling 
laws of SU average capacity with respect to F^ and F t f . ■ 

Further, it can be also shown that for the multiple PUs 
case, the average capacity of the mth SU is converging to 
the lower bound on average capacity for the single PU case as 
N, F — > oo. Assume without loss of generality that an infinite 
number of subcarriers F is available. Because the orthogonal- 
ity of PUs' subcarriers is maintained, then 2~2 n =i ~ F as 
F, N — > oo. Hence, 



lim E \CL 

N,F-hx L 



lim 



F S 

m 

i~F~ 



' N 

E 

.n=l 



F, P E 



C 



I.n 



-FfE [C, 



NI] 
rn.i] 



= F S E 



C 



I.n 



where it is assumed that 
power. Thus, E C 



I.n 



all the PUs have the same transmit 
is the same for all N PUs. 



6 



IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. XX, NO. XX, XXX 2012 



Corollary 4: The average capacity of the ?7ith secondary 
user in the presence of a single PU converges logarithmically 
to F^E [C^A as F increases towards infinity: 



> 

with log(F) 



F* E [C%] 



(3) 



Proof: The proof is given in Appendix [C] ■ 
Using similar steps, one can readily obtain the bounds and 
the scaling laws of the SU average capacity in the presence 
of multiple (N) PUs in the primary network. 



B. SU Capacity Analysis over Rayleigh Channel Fading 

In this section, the SU capacity over a Rayleigh channel 
fading model is investigated. Thus far, the CR capacity studies 
in the literature have mostly assumed two types of PUs' 
interference constraints on the SU transmit power: the peak 
power interference constraint and the average interference 
constraint J6}, (7). The peak power interference constraint is 
adapted in this work, and an adaptive scheme is used to adjust 
the transmit power of SU to maintain the QoS of PUs. Hence, 
the transmit power of the mth SU corresponding to the ith 
subcarrier is given bjj^J 



jVj 1 J Pm : i ) ^ F m ih m p^i 

m ' 1 = I ^ * <P h 
> ;> 

itl\i\< fVi, -. 



Let X m ,i — h>m,iPm ^ tnen tne received SINR of the mth 



for i = 1,...,F 
Let A mj i = h. 
SU's ith subcarrier is 



7^-, forn=l,...,N, (4) 



where I^i = Pn,i9ns,i stands for the mutual interference 
caused by nth PU on the ith subcarrier. In Q, S^™ represents 
the SINR in case when subcarrier collision occurs. Therefore, 
when there is no collision, i.e., the subcarrier is not being used 
by two users, there is no interference caused by PUs. Hence, 

The CDF of A m , can be obtained as follows f30): 



F x .{x) = F hmp . 



p 



Fojx) 



F 



& 3 \hr, 



p 



where ■d 1 = h m ^P m<l and d 2 = ^ih mji /h mPfi , with their 
corresponding PDFs given by f^(x) = e^ x ^ i and 
fy 2 (x) = ^i/(x + *i) 2 , respectively. Hence the CDF and 

5 Notice that due to the random allocation, the SU transmit power is adapted 
(regulated) considering the worst case scenario, as if all the subcarriers in 
the primary network are utilized by PUs. This condition assures the QoS 
requirements of PUs. 



the PDF can be expressed, respectively, as 



1-e 



= 1-e" 
dF Xmi (x) 



e -^W (5) 



dx 



P, 



1 - e 



.X 2 + %X - * 4 P m>i 



(*i +X) 2 



(6) 

Similarly, by using a transformation of RVs, the PDF of 
^mi w i tri fi p (v) — e~ y / PnA /P n ,i can be expressed as ||7j 

F s i,n (x) = Pr (A m>i < x + rj)) 



F\ m < ( x (v + r n))fi p Xy)&y- 



(7) 



Plugging (|5]l into (|7j, it follows that 



jcPl i ~ F, 



-L 



1 



x r 0, r? + 



p m .i 



xPn 

X 



pp. 



where the upper incomplete Gamma function is defined as 
r(x, y) — J t x ~ 1 e~ t dt, and the derivation of CDF yields 
the PDF V 

XT)Pn,i + P m ,i{V + Pn,i) ( 3- ,\ -25±^i 



(xP n ^i -\- Pm,i) 



x e 



(*< + xP n>i )T 0, jj 



11 e p "..« 
*,\ / 1 



x 3 p2 



xP n 4x 2 r]P n ^ - ViP m ,i) c -( v+ ^)( 7s 2- + 1 ^-) 

{xr) + ^i){xP nA + P mti ) e 

(8) 

Similarly, when there is no primary interference using |6]) 
and the transformation f S Ni(x) — f]f\ m ^x), it follows that 



fs™ ( x ) 
p '. 



1-e 



^ / ( T yx) 2 + ^??x-^P m , j 
V (*i + ?yx) 2 



(9) 



7/J' 



17 3; + *^ 



(10) 



and the CDF is given by 

F S M I (x) = l-e-^ 

"*•* Wj + rjx 

Finally, the desired expressions for the PDFs of C^™ and 
can be obtained by transforming the RVs as follows: 



f c i,n (x) 



dy 
dx 



= e»/ s x.»(e»-l) ) 
„=e.-l m - (11) 

/c"'» =e x f s *i. (e x -l). 
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Using Definition |4j for any arbitrary mth SU and multiple 
(N) interfering PUs, the instantaneous SU capacity with 
subcarrier collisions is given by 

C m = £log (l + +■■•+£ log (l + S^ 



i=l 



i=l 



for n = L 2, . . . , N, and var(x) denotes the variance of x. 
From |7|, using ([8j)-(fTTj), the average capacity of C^™ and 
can be expressed, respectively, as 



E 



C 



I.n 



xf c i,n (x)dx 



5>g(l+<i) 




(12) 



N k nm fe/m 

/ j m,i ' / ^ m,i ' 
n— 1 z— 1 i—l 



CI 



There are two types of well known methods available to 
evaluate the distribution for sum of variates, namely, the 
characteristic function (CF) and the moment generating func- 



tion (MGF) based methods [23]. Unfortunately, by employing 
these methods, it is often hard and intractable to obtain 
explicit closed form expressions for the PDF and CDF of 
SU capacity in ( [T2] > from (|8]>-([TT|. Even if we obtain, it 
will hardly provide any insights because of the complicated 
expressions. Therefore, in order to sum up the rates for the 
cases of interference and no-interference, we will approximate 
the PDFs of Cjfi and C%\ using a Gamma distribution. There 
are important properties of the Gamma distribution that are 
suitable for approximating the PDFs of the variables C^ 7 ] 
and C^\. First, the sum of Gamma distributed RVs with 
the same scale parameters is another Gamma distributed RVs. 
Second, the skewness and tail of distribution are similar for 
the whole range of interest and are determined by mean and 
variance (26). Last but not least, Gamma distribution is a 
Type-Ill Pearson distribution which is widely used in fitting 



positive RVs [24|-|26|. In addition, since Gamma distribution 
is uniquely determined by its mean and variance, we employed 
the moment matching method to the first two moments: mean 
and variance. 

Definition 5: X follows a Gamma distribution, X ~ 
G(a, /3), if the corresponding PDF of X with scale and shape 
parameters, /3 > and a > 0, respectively, is given by 



fx(x) 



exp 



(->) 



U(x), 



f3 a T{a) 

where U(-) denotes the unit step function, and the Gamma 
function is defined as T(x) = J Q t x ~ 1 e~ t dt. 

Since the mean and variance of Gamma distribution are a/3 
and a/3 2 , respectively, mapping the first two moments with 

yields 



the PDFs of and C. 



NI 



NI 



var 



a 



I.n 



var 



I,n 



E 



a 



I,n 



l0g(l + X)fgl,-n (x)dX 



1 - e 



1 - 



x r o, 



p 

1 m.i 



J n..r 



P 

1 n, 



r o, u 



-r(o,-^ ie- 



x(l + x) 



-dx, 



(13) 



and 



EK'J = / xf C NX i (x)dx= / ]og{l+x)f s m(x)dx 



r o, 




(14) 



— — r o, — 

V P, 



The variance of C^™ is given by 



C' 



E 



(^") 2 1-(e [eg; 



where the second moment of C^j" is expressed as 



E 



= / [\og(l + xff sI , n (x)da 



21og(l + at) 



1 + x 



1-iWz) 



o 

N p 
3 = 1 



21og(l + s J ) 

' 1 + Sj 



l-F s r, n ( Sj ) 



where the second equality is obtained by using integration by 
parts [|7|. The resulting integral is readily estimated by em- 
ploying Gauss-Chebyshev quadrature (GCQ) formula, where 
the weights (wj) and abscissas (sj) are defined in [31, Eqs. 
(22) and (23)], respectively. The truncation index N p could be 
chosen to make the approximation error negligibly small such 
as N p — 50 for a sufficiently accurate result. 
Similarly, the variance of is expressed as 



Nil 



E 



NI\ 
m.i ) 
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where the second moment of is calculated as follows 



E 



i=i 



21og(l + Si ) 

1 + Sj 



Therefore, using the Gamma approximation, the capaci- 
ties are approximated as C^;" — S^iAi) an d ~ 

In Figure [2] the exact and approximative expressions of 
f c i, n (x) and f C Ni (x), including the simulations results, for 
different system parameters are shown. It can be observed that 
the approximation is very close to the exact results. 

Since both and are i.i.d. for given k nm , the con- 



ditional characteristic functions for the rate sums 12* 
and X^=7 can t> e expressed as follows 



fen 



m r l 
1 ^m,i 



= (1 - J^ NI 



-Ceifen 



$ c «/(o;|A;„ m ) = ^$ c nj(w) 

where ^^.n (a>|fc nm ) and $ c n; (w|fc nm ) are the characteristic 
functions of f c i,n(x\k nm ) and / c i»j (x\k nm ), respectively. 
Using the nice feature of the Gamma distribution that the sum 
of i.i.d. Gamma distributed RVs, with the same scale param- 
eters (/3) is another Gamma distributed RV, the conditional 
PDFs are expressed as follows 



fck n \k n , 



(x\k n 
(x\k n 



t m k< <^ 



(15) 



In ([12), even though the conditional PDFs of C^ n and C^ 1 
are obtained, to find the PDF expression for C m , we first need 
to evaluate the PDF of Cf n , and then the PDF of its sum 
with C^ 1 . At this point, one needs to be aware that there are 
N + 1 terms in ( [12) , and each follows a Gamma distribution 
where the shape (a) and scale (fi) parameters can be arbitrary. 



Therefore, the aforementioned feature of Gamma distribution 
for a sum of Gamma variates cannot be employed here. 

Expressions for the PDF of sum of Gamma RVs are derived 
by Moschopoulos J27) , Mathai ]32) , and Sim J33] , In addition, 
constraining the shape parameters to take integer values^] and 
be all distinct, by using the convolution of PDFs Coelho p4| 
and Karagiannidis et al. ]35) , or partial-fractions methods 
Mathai p2[ , derived an expression for the PDF of a sum 
of Gamma RVs. Nevertheless, Moschopoulos PDF provides 
a mathematically tractable solution that it does not restrict the 
scale and shape parameters to be necessarily integer-valued or 
all distinct (32). Therefore, the following theorem will help us 
in this regard. 

Theorem 2 (Moschopoulos, 1985 j^j): Let {X s }f =1 be 
independent but not necessarily identically distributed Gamma 
variates with parameters a s and f3 s , respectively, then the PDF 
of Y = X) s =i X s can be expressed as 



fy(y) = LI 



k=Q 



5 

12 a s 

s=l 



(16) 

where /3i = min s {/3 s }, and the coefficients Sk can be obtained 
recursively by the formula 



5 = 1 



-, fe+i 

— y 

i— i 



Jfe+l-i 



for k = 0,1,2,... 

Proof: See (27). ■ 
The Moschopoulos PDF provides a nice and tractable rep- 
resentation of sum of Gamma variates in terms of a single 

6 If the shape parameter is an integer, Gamma distribution is referred to as 
Erlang distribution. 
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(a) 



(b) 



Fig. 3. Moschopoulos (a) PDF j 1 6| and (b) CDF j 1 8| for ?t = 25, and the total number of Gamma distributed RVs in sum as S = 4 and 5 = 2. 



Gamma series with a simple recursive formula to calculate the 
coefficients. This representation is applicable for any arbitrary 
shape parameters {a s }f =1 and scale parameters {/3 s }f = i 
including the possibility of having some of the parameters 
identical. 

The CDF of Y can be obtained from the PDF as Fy (y) = 
fy{x)Ax. Therefore, 



Mv) = 11 



s=l 



a s oo 

E 



3fe 



exp 



x 

Ti 



s 



dx. 



(17) 



The interchange of summation and integration above is 
justified using the uniform convergence of (16) (see e.g., pj\ 
for a rigorous proof). From (36), we can simplify (17} by 
using J Q x"~ 1 e~> J,x dx = \i~ u ^{v,\xu) for !ft > 0] |37] pg. 
346, Sec. 3.381, Eq. 1], where 7(-, •) is the lower incomplete 
Gamma function and is defined as j(x,y) — Jq t x ~ 
Hence, 



Hi. 



s 

Mv)=W 

s=l 



n 



7 



s 

s=l 



fc=0 



fc=0 




(18) 



where V(-,-) is the regularized (also termed normalized) in- 



complete Gamma function and defined aa^T^ 



1 - 



r(a,z) 
T(a) 



T(a) 



For practical purposes, based on the required 



7 For integer values of 5Z?=i a s + k, using 1371 Eq. 8.353.6] the 
regularized incomplete Gamma function can be further simplified to 



■PITT. 



l-exp(-y/&)£tS 



Ef=i"s+fc-i i 



accuracy of application one may use the first /i, i.e., k = h — 1, 
terms in the sum series ( [15} , The expression for truncation 
error is given in p7) . In Figure [3] the Moschopoulos PDF 
and CDF are shown for S = 4 and 5 = 2 where only 
the first 25 terms in the infinite sum series, i.e., h = 25, 
are considered. One can observe that the Moschopoulos PDF 
and CDF perfectly agree with the simulation results for same 
values of a and f3. Since in our system model, with some 
probability the transmit power of PUs P„ ,i for n = 1, . . . , N, 
can be the same, which means that the corresponding of n and 
/3 T { are the same. Such a scenario can arise when the PUs are 
at the same distance from their corresponding common PBS. 

Recall that from ( p~2| > and (IS) , we have to evaluate 
the PDF of the sum C^ 1 + C$ + ■ ■ ■ + C^ N + C^ 1 , 
for a given number of set of subcarrier collisions k m = 
[kim, k 2m , ■ ■ ■ , k Nrn , kf m ]. Recall also that C/„ and C^ 1 
are Gamma distributed and independent but not necessarily 
identical. Therefore, the conditional PDF of their sum can be 
expressed by means of Theorem |2] as given in (19) , where 
/?min = min{/3( , /3| , . . . , j3L, (3 N *}, and the coefficients 8k 
are obtained recursively as follows: 

<5o = l 



S k = 



1 



fc+i 



, NI 



N 

E< 

i - 



1 



P 



NI 



Sk+i-i for k = 0,1,2,.. 



Now, the PDF of C m can be found by averaging over the 
PMF of subcarrier collisions as follows: 

l {x, k T , 



fc m (x)=J2f< 



C m ,K„ 







E 



fc m |K m (^|k m 



(20) 



Plugging ([T| and (jT9j into (20) results the PDF, and is given 
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fc m \K m [x\K*) = [ im ) 11 hr Z^ — — - - — - ( n y (19) 



n=l 



fc.w=EE-EE t ' u s n w n 



Pmin 




fcjVm & /' 



(21) 



The outage probability is a common performance metric 
in fading environments. Hence, here we consider the outage 
probability of SU capacity in terms of the following measure: 

Pc2(<Pth.) =Pr(C m <ip th ) 

<Ptb 

fc m {x)dx, 



which is the CDF of the SU capacity over the outage threshold 
</?th [dB], Using {L8]i and the CDF of C m can be 

expressed as 




^ SkV ^ a^knm + a NI k fm + k, 



k=0 



V. Asymptotic Analysis of Multiuser Diversity 

In this section, the gain of multiuser diversity by employing 
opportunistic scheduling is investigated. In conventional sys- 
tems, the multiuser diversity gain is attributed to channel gains 
only. However, in the proposed scheme, we additionally benefit 
from the randomness of the number of subcarrier collisions. 
Assuming all M SUs are accessing the F available subcarriers 
to randomly allocate their subcarriersj^jthe SU, which provides 
the best instantaneous capacity, is selected as: 



C n 



max Cn 
me[i,&f] 



For fairness in the selection phase of the best SU, assume 
that each SU's data rate is the same, i.e., each SU requests for 
the same number of subcarriers, — F s , m = 1, . . . ,M. 

8 It is assumed that no collisions occur among the subcarriers of SUs. 



Then, by using order statistics, the PDF of C max is expressed 



fc BM (x) = Mf Cm (x)F Cm (x) 



M-l 



(23) 



Plugging ([21} and ([23) into g3), the PDF of C max can be 
obtained. Nonetheless, using /_ x fc mm {x)dx is intractable 
to find the mean of C max . Even if we can carry out such 
a calculation, it will hardly provide any insights to fully 
understand the impacts of the main parameters on the capacity 
using the resulted expression. Therefore, we asymptotically 
analyze the capacity to understand the effects of system 
parameters and multiuser diversity gain in CR systems with 
spectrum sharing feature. 

Theorem 3: As the number of SUs M goes to infinity, the 
average capacity of C max converges to 



where E\ = 0.5772... is Euler's constant |38|, and aj\/ = 
[MfcJbAiT 1 and b M = F^(l - 1/M). 

Without loss of generality assuming a single PU case, i.e., 
n E [1, N], then b^j is given by 



F' 1 



1 



1 

M 



Qj2 SkV ^ ( A + fc ' 



k=0 



J_ 



where "P" 1 ^, •) stands for the inverse regularized incomplete 
Gamma function. Unfortunately, there is no closed form 
expression for this special function. Therefore, it can be 
evaluated numerically by using build-in functions in some 
well-known computational softwares such as MATLAB® and 
MATHEMATICA®0 Additionally, /3 min = mm{l3^/3 NI }, 

9 It is worth to note that by using [39 6.5.12 & 13.5.5], the regu- 
larized incomplete Gamma function can be approximated as V (u, v) = 
> iFi <u:l + u; —v) = Sr. , as v — > 0, where ifi (•;•;•) is 

ur(v) J- V * ) > u Tt u) \ ' ' / 

confluent hypergeometric function 1 40 1. Hence, its inverse can be obtained. 
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A = alk r , 



t NI kf m and Q takes the form: 



Q = 



F 
F s 



E 

fc„ m =0 



± n 
h 

n i i.i >i 



F-F p 



rmin 



where in considering a practical scenario, it is assumed that 
F^ + F p < F and F^ < F% ■ Hence, the support region for 
the number of subcarrier collisions is k nm = 0, 1, . . . , F^. 
Proof: We start with the following Lemma. 
Lemma 1 (Distribution of Extremes ^38^ ): Let z±, . . . , Zm 
be i.i.d. RVs with absolutely continuous common CDF, F(z), 
and PDF, f(z), satisfying these conditions: F(z) is less than 1 
for all z, f(z) > and is differentiable. If the growth function 
g(z) = (1 — F(z))l f(z) satisfies the von Mises' sufficient 
condition: 

lim g(z) = c> 0, (24) 

z— >oo 

then F(z) belongs to the domain of attraction of the Gumbel 
distribution. In other words, [maxi<fc<M z m — &m]/om con- 
verges in distribution to the Gumbel-type limiting distribution: 



G(x) = cxp (— t 



— oo < z < oo 



Thus, the maximum of M such i.i.d. RVs grows like &m, 
also termed position parameter. The parameter bj^ is given by 
i>M = — 1/M), and the scaling factor au is given by 

a M = g(b M ) = [Mf(b M )r 1 - 

The PDF and CDF of C m for a single PU are given, 
respectively, by 



°° Si.x A+k - 1 e- x /^ 



F c ijx) = Qj2s k V[A 

k=0 



k. 



Pn 



(25) 



(26) 



with the coefficients calculated iteratively as 

6 = l 

k+l 



— E 



+ a NI k fm 



a k 



1 



A, 



Bn 



8 



AT J 



4+i-,: forfe = 0,l,2,. 



From Lemma [T] plugging ((25) and |26]) into |24} yields 

1 - Fm (x 



lim 

x— f oo 



A, 



> o. 



(27) 



The respective intermediate steps in the evaluation of ( |27| ) 
are depicted in Appendix [D] Hence, it belongs to an attraction 
domain of Gumbel-type with limiting CDF: 



F c 



. (x) = exp I— exp 



x - b M 



Then, the limiting PDF of C max is 



1 f x-b M \ ( ( x - b M 

— exp exp — exp 

a-M V au J \ \ a M 



Therefore, using E [C max ] = xfc lnllx {x)dx, the desired 
result can be readily obtained. ■ 

In the proof stage, it came to our attention that, to the 
best of the authors' knowledge, there is no result reported in 
the literature for the limiting distribution of RVs that follows 
Moschopoulos PDF. Therefore, the following novel result can 
be stated. 

Corollary 5: Let {X r }^ =1 be the set of 1Z i.i.d. RVs 
that follow Moschopoulos PDF and CDF (27), and Y = 
max {Xi, X2, ■ ■ ■ , X-ji}, then the limiting distribution of the 
CDF of Y belongs to the domain of attraction of Gumbel 
distribution as 1Z converges to infinity. 

Proof: It is immediate to see this result from the results 
presented in the proof of Theorem^ ■ 

The results obtained so far will help us to asymptotically 
analyze the scheduling of SUs' subcarriers in the following 
section. 



VI. Centralized Sequential and Random 
Subcarrier Allocation 

A. Sum Capacity of SUs with Multiuser Diversity 

In this section, a cognitive communication set-up involving 
multiple SUs that assume a random allocation method is 
studied. Recall that due to random allocation scheme, there 
can be the collisions among the subcarriers of SUs in ad- 
dition to those that are used by PUs. These collisions will 
decrease the system performance severely. To overcome this 
challenge, we propose an efficient algorithm that sequentially 
and randomly allocates SUs' subcarrier sets in a centralized 
manner by maintaining the orthogonality among the allocated 
subcarrier sets. Such an assignment can be thought of as 
the downlink scenario where the SBS performs the random 
assignment of subcarriers. Furthermore, to benefit from the 
multiuser diversity gain, the opportunistic scheduling method 
is employed in the algorithm [See Table [I], where it is assumed 
only a single PU. The multiple PUs case is a straightforward 
extension. In the selection step of the best SU, to preserve 
the fairness among the users, it is assumed that the data 
rate requirements of all SUs are the same, i.e., the individual 
numbers of subcarrier requirements are equal. 

The algorithm can be summarized as follows. A randomly 
chosen set of subcarriers Ff from the set F is assigned to the 
available SUs. The first SU, which provides the best capacity, 
is selected among M SUs, then the selected subcarriers F^ 
(total of collided and collision-free subcarriers) are removed 
from the set F. In the next stage, another randomly chosen set 
of subcarriers Ff from the updated set F is allocated to the 
rest of SUs. The second best SU is selected among the M — 1 
SUs, and similarly the subcarrier set F is updated by removing 
the new set F^-. This sequential selection continues until it 
reaches the total number of the best M SUs, with M < M. It 
is evident to observe that the multiuser diversity is attributed 
only to the randomness of the channel gains. Furthermore, 
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TABLE I 

ALGORITHM: CENTRALIZED SEQUENTIAL AND RANDOM SUBCARRIER ALLOCATION 



1) Initialization 

• Assume = F s Vm £ [1, M] and a single PU is available, n = 1. 

• Set the number of available subcarriers to F and index t = 1. 

2) Subcarrier assignment step 

• Randomly sample a set of subcarriers, F R , with cardinality of F s from set F: k nm ~ HYPG(F S , F,f , F). 
. Assign the set F t R to all M - t + 1 SUs. 

3) Capacity calculation step 

• For m = l,...,M-i + l, SUs evaluate their capacities with the given random set of subcarriers: C m F t . 

• SUs send feedback for the calculated capacities to the central control entity (SBS or CR Network Manager). 

4) Selection step 

• Choose the SU that provides the best capacity: 
lit =1 then ml = argmax (CmjF/ 1 ) 

m£[l,Af] 

else m* — argmax (C m \F R ) for t = 2, . . . , M. 

m£[l,M] \ [m£ i m t — l] 

5) Updating the subcarrier sets step 

• Remove the sampled (total of collided and collision-free) subcarriers from the available set of subcarriers: 
F -S- F- F t R . 

. Set t <- t + 1 and go to Step 2 until t = M. 

6) 5wm capacity evaluation step 

M 

• Compute sum capacity of SUs: C sum = XI ^m* • 

£=1 



some of the essential points in the algorithm can be highlighted 
as follows: 

• In step 2: The PMF of the number of subcarrier collisions 
follows a hypergeometric distribution due to the random 
selection of subcarriers set from the available set of 
subcarriers F. 

• In step 4: The selection of the best SU is performed based 
on the capacity feedbacks obtained from the SUs. 

• In step 5: Removing the randomly sampled subcarriers 
F R from the available set of subcarriers F means that 
both collided and collision-free subcarriers are subtracted 
from set F (since Ff — k nm + kf m ), i.e., F F — 
F t R ^ F <- F- k nm - k fm and F£ <- F£ - k nm . In 
other words, since F = Fjf + Ff, where Ft stands for 
the number of free subcarriers, the subcarriers that are 
occupied by the PU Fjf in the set F are automatically 
updated when the randomly sampled set of subcarriers 
F 4 is removed from the set F. Hence, the orthogonality 
among the subcarriers of SUs is maintained. 

Theorem 4: The sum capacity of M selected SUs in the 
centralized sequential and random scheduling algorithm for 
M ^> M is approximated 10 ! by 



E[C a , 



ME [C„ 



and as M — > oo, it converges to 

E [C sum ] = M [b' M + E ia ' M ] , 
where ml is the index of the first selected best SU and defined 

10 Since E \c m *\ > E \c m *] ,Vj £ \l,M], it can also be considered 
as a tight upper bound for M > M as E [C sum ] < ME \ C„ 



as ml = argmax C m . Further, a' M and b' M can be readily 

m£[l,M] 

obtained by following the same approach as in Theorem [3] 
considering the fact that the multiuser diversity is only ascribed 
to channel randomness not the random subcarrier assignment. 
It is noteworthy to state that the sum capacity scales linearly 
with the number of selected SUs. 

Proof: The scheduler selects the SUs according to the 
following rule: 

to* = arg max C m for j = 2, . . . , M , 

mtE [1,M]\ [m* j 

which means that the selected SU(s) are ignored in the 
selection step of remaining users. Then, the sum capacity of 
selected SUs is defined as 

M 

Csum = ^ ] ' (28) 

For large M such that M ^> M, it is immediate to observe 
that 



E 



1, M 



(29) 



This approximation is valid since removing the selected 
SUs does not considerably impact the mean of the rest of the 
selected SUs for M ^> M, i.e., the maxima of M RVs and 
M—M RVs are approximately the same for M 3> M, so their 
averages are approximately the same. Hence, plugging p9] ) 
into ( |28| i yields the desired result. ■ 

B. Sum Capacity of SUs without Opportunistic Scheduling 

In order to investigate the performance of our proposed 
algorithm due to multiuser diversity gain, the performance 
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Fig. 4. SU mean capacity versus the transmit power P m ,i with different IT 
*i values for = 20, F 7 f = 30, F = 128 and P n>i = 10 dB. 



Fig. 5. SU mean capacity versus the interference temperature >T/; with 
different transmit power P m>i values for F^ = 20, F% = 30, F = 128 
and P n> , = 10 dB. 



of the centralized sequential subcarrier scheduling without 
employing the opportunistic scheduling method is analyzed 
in this section, i.e., the multiuser diversity of SUs is not 
maintained. Therefore, the sum capacity of any arbitrarily M 
selected SUs (among M SUs) can be expressed as 



M 

E 



Recalling the upper and lower bounds on the average 
capacity of a single SU C m , one can conclude that the average 
sum rate of the SUs scales linearly with the number of selected 
SUs (M). Mathematically speaking, E [C s a um ] = ME [C m ]. 

During the sequential scheduling of SUs' subcarriers, the 
PMF of the number of subcarrier collisions can be obtained 
as a special case of the following result. 

Proposition 2: The PMF of the number of subcarrier colli- 
sions for the mth SU in the presence of N PUs, when assign- 
ing the subcarriers sequentially to preserve the orthogonality 
between SUs' subcarriers, is given by 



EE'" E P( k l,k 2 ,...,k m ), 



where the joint PMF is 



p(ki,k 2 

Ff 
k fl 



F 
if 



! k m ) — 

-1 N 



n 



£• 



FS 



• »i 

-l 




N 

n 



F 1 



h 

n n r 



The mean and support of k nm are given, respectively, by 

/ m-1 \ 

[Fn- E Mknj] 

E [& nm ] = 



3=1 



m— 1 

F- E Ff 

3 = 1 



N 



m—1 



kf m = F,i and k nm £ ( F? n + F n 



- k nj -F\ ,...,mm I if ,F? - ^ k n 

3 = 1 / [ 3=1 

Proof: The proof is given in Appendix [E] ■ 

VII. Numerical Results and Simulations 

In this section, numerical and simulation results are pre- 
sented to confirm the analytical results and investigate the 
impact of various system parameters in CR spectrum sharing 
networks. First, the effect of peak transmit power of SU P m ,i 
(in dB) on the average capacity (in nats per second per hertz) 
is shown for different values of IT values "J/j in Figure |4] 
Unlike the conventional systems, the SU average capacity is 
here saturated after a certain value of peak SU transmit power 
because of the IT constraint in spectrum sharing systems. In 
Figure [5] the SU mean capacity against the IT constraint is 
presented. It turns out that the analytical results agree well 
with the simulation results. The results shown in Figures [4] 
and [5] are in the presence of a single PU, i.e., n 6 [1,-^V], 
and the number of subcarriers in sets F, F 7 ^ and F*~ are 
chosen arbitrarily^] A common observation for both Figures [4] 
and [5] is that the saturation level of capacity increases as the 
IT constraint relaxes, and the capacity keeps growing until 
a saturation point as the transmit power of SU increases as 
expected. It can also be underlined from Figure [4] that the 
capacity gain due to relaxation in the IT constraint disappears 
at low SU transmit power. Therefore, in the high transmit 
power or SINR regime, the impact of IT relaxation differs 
significantly. Similarly, the same effect can be observed for 
the results in Figure [5] 

"The unit AWGN noise variance is used (77 = 1) in all the following 
figures. 
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SU mean capacity versus the transmit power P m ,i with different 
of PUs N for F^ = 20, F = 128, F% = 10, P n>i = 5 dB for 
...N and *i = -5 dB. 



Consider now the practical scenario when there are multiple 
PUs available. Therefore, the number of free subcarriers in 
the available set F is smaller than that of the single PU 
case. The SU mean capacity against peak transmit power 
P m ^i in the presence of multiple PUs is shown in Figure [6] 
In order to illustrate the effects of multiple PUs, during the 
simulations, it is assumed that the number of subcarriers and 
the transmit power of all PUs are the same, P n i = 5dB 
and F 7 f = 10 for n = 1,...,N, respectively. Since the 
number of subcarrier collisions in the presence of multiple 
PUs follows a multivariate hypergeometric distribution, the 
multivariate hypergeometric random variates are generated by 
using the sequential method given in pT[ p. 206]. It can 
be observed that increasing the number of PUs degrades the 
performance of SU as expected. In addition, as the number 
of PUs decreases, i.e., the number of unoccupied subcarriers 
increases, the average capacity of SUs converges to the upper 
bound, where all SU's subcarriers are collision-free. On the 
other hand, the lower bound of the average capacity indicates 
that all SU's subcarriers are colliding. 

Figure [7] shows how the SU average capacity scales with the 
number of subcarriers in sets F and F£, respectively, where 
the single PU case is assumed. As the number of available 
subcarriers increases for a fixed number of SU's and PU's 
subcarriers, the SU mean capacity asymptotically converges to 
the limit point given in Q, where the rate of convergence is 
logarithmic. It is immediate to see that the SU average capacity 
scales as 6 (1 + 1/F), 6 (F%) and 6 (l - F£), as proved 
in Corollary^ 

The performance of the proposed centralized algorithm 
with and without the opportunistic scheduling is simulated 
and shown in Figure [8] The results for the algorithm with 
multiuser diversity are in the presence of M = 10 and 
M = 40 SUs and among them M = 5 SUs are selected 
using the opportunistic selection method. Also, M — 5 
SUs are selected when no opportunistic scheduling method 




20 40 60 



80 100 120 140 160 180 200 

Number of Subcarriers 



Fig. 7. SU mean capacity versus the number of subcarriers F and F^ , for 
P n>i = 5 dB, P mJ = 10 dB and = -5 dB. 
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— B — Algorithm [M - 10] 
— 6 — W/o opportunistic scheduling [M = 1 0] 

W/ collisions between SUs [M = 1 0] 

W/ collisions between SUs [M - 40] 
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Fig. 



Sum capacity of M = 5 selected SUs versus the transmit powers 
10, m = 1, . . . , M, F£ = 40, F = 100, P n>i = 10 dB 
and *, = dB. 



is employed. Note also that without opportunistic scheduling, 
the number of SUs M does not affect the sum capacity of 
M selected SUs. Therefore, this scheme is not plotted for 
different numbers of SUs (M). One can observe that the effect 
of multiuser diversity manifests into the fact that an increase 
in the number of SUs M results in higher capacity in the 
proposed algorithm. Furthermore, in order to reveal the impact 
of collisions between SUs subcarriers on the sum capacity of 
any arbitrarily M selected SUs, we simulate the performance 
of M = 5 selected SUs in the presence of M = 10 and 
M = 40 SUs in the secondary network when no centralized 
algorithm with opportunistic scheduling is employed. In other 
words, the orthogonality among the subcarriers of SUs is not 
maintained, and the multiuser diversity gain is not exploited. 
Hence, there can be the collisions between the subcarriers of 
any SU with the rest of the SUs in the secondary network in 
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addition to those that are utilized by PU. This scheme could 
be considered as the worst case scenario, where the collisions 
among the SUs' subcarriers severely affect the performance 
due to high probability of interference level among SUs as 
shown in Figure [8] 

VIII. Conclusions 

This paper studied the performance of OFDM-based CR 
systems with spectrum sharing feature using a random sub- 
carrier allocation method. The subcarrier collision models for 
single and multiple PU(s) are shown to assume univariate 
and multivariate hypergeometric distributions, respectively. 
The expressions of SU average capacity for both general 
and Rayleigh channel fading models are presented. It turns 
out that the closed-form expression for the instantaneous 
SU capacity in the presence of Rayleigh channel fading is 
intractable. Therefore, the Gamma approximation of the SU 
capacity expression is obtained by employing the moment 
matching method and Moschopoulos PDF representation for 
a sum of independent but not necessarily Gamma distributed 
RVs. Through the asymptotic analysis of SU mean capacity, 
it is found that the capacity scales with the number of 
subcarriers as 0(1 + 1/F), (F%) and G(l-F£). The 
asymptotic analysis of capacity assuming an opportunistic 
selection method is investigated by using extreme value theory. 
When multiple SUs are randomly allocated the subcarriers, the 
primary issue that causes drastic performance degradation is 
the collision(s) among their subcarrier sets. In order to prevent 
such a situation, a centralized algorithm was developed to 
sequentially assign orthogonal subcarrier sets to SUs based 
on a random allocation scheme while benefiting from the 
multiuser diversity gain for maximum SUs sum rate. Besides, 
it is found that the extreme value limiting distribution of RVs 
that follow the Moschopoulos PDF belongs to the domain of 
attraction of the Gumbel distribution. 

Appendix A 
Proof of Theorem[T] 

According to Definition |4] to evaluate the average of sum 
capacity of the SU with subcarrier collisions, we have to 
average a random sum of RVs with the set of i.i.d. RVs C^ 7 ] 
and as follows: 



E [C* ] = E 



E 



■E 



E n I,n , \ n NI 

i 

k fm 



E 



K — h 



iNI 



= E 



j=l 

E [k nm E [C 1 ^ 



■ E 



+ E[fc /m E [C$]] 



where the rule of iterated expectations |42| p. 55, Theorem 
3.24] also known as tower rule, E [X] = E [E [X\Y]], is ap- 
plied, and the conditional expectations with respect to k nm ~ 

HYPG(F^, Ff, F) and k fm ~ HYPG(F^, F — F% , F) are 
used. 

Furthermore, k nm and C^' n 



are independent, and so are 



and C^\. Then, we have 



E [C* ] = E [k nm ] E 



C 



■E[k fm ]E[CZ%. 



It is also worth to note the relation between the two sums 
in the first equality that they are independent conditioned with 
the given k nm and hfm (since kfm — — k nm ). Taking into 
account the means of k nm for n£ [1,N] and fc/ m , it follows 
that E [k nm ] = F^FP/F and E [k fm ] = Fg{F - F£)/F, 
which yield the desired result. 

Appendix B 
Proof of Corollary[T] 

Following the same approach as in Appendix [A] the average 
capacity in presence of N PUs can be obtained by using ( p"2] ) 
and the properties of multivariate hypergeometric distribution 
given in ([1} with the means of k nm and kf m expressed 

as E[k nm ] = FgF r f/F, n = 1,...,N and E[k fm ] = 



Let xi 



F„ E 



= E ^m 

can show that 



Appendix C 
Proof of Corollary|4] 

F s 

-frXi 



E [C%<] 



X2 



e [a 



FmX2- Using Definition 



NI 

m,i 



i 



and 



lira 

F->oc 



lim 

F-KX) 



F S 



F» 

F+lXl 



F s 



X2 



F+lXl 



F^ 



X2 



-frXi 



F^ 



X.2 



= 1, 



and 



lim 



°m,F+l 



Ff n X2 



°m,F 



F s 



X2 



lim 

F->oo 



F+l 



CI + F t n 



X2 



F S 



X2 



FS 

-#Xi 



F s 



X2 



F s 



X2 



Hence, 



as F 



is logarithmically convergent to F^E [C^a] 



Appendix D 
Evaluation of Limit in Equation ([27]) 

In the evaluation steps, given in pO) , since fai (x) — > and 
Fqi (x) — > 1 as x — > co, first, L'Hopital's rule is applied, and 
then due to the uniform convergence and the positive terms, 
the interchange of limit and infinite sum is viable. Lastly, 
because the resulting expression is of polynomial type, only 
the highest-order terms are considered. 
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Iim — — -t — = lim 

x— >oo Jqi [x) x— foo 



OO t 
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k=0 



1 r ( A+fc '^) ' 

1 r(A+fe) 



fc=0 



/ 9^+ /e r(A+fc) 



2 E Ail 



<5 fc a; A + fc - 1 e- x /' 3 min 



= G r(A+fc) 



C/(x) 
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lim 
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k=0 



r(A+fe) 
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2 E fl >- 
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lim 

Ik— too 
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Ptt' k r(A+u) 
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(30) 



Appendix E 
Proof of Proposition|2] 

Recall from ([]} that the PMF of first SU is given by 



Z)(5)~!r 



Assuming the orthogonality between subcarriers, given k x 
the conditional PMF of second SU is a multivariate hyperge- 
ometric distribution, is described by 

p(k2|ki)= l J l f§ ) n 



p - k i 

2 / "J \ fc «2 
AT 



where 1 T = [1, 1, . . . , 1] T e Z N+1 and l T ki = E„ =1 *Vu + 
fc/i = Ff. Similarly, for the third SU the conditional PMF 
for the number of subcarrier collisions is 



p (k 3 |ki,k 2 



F f - k fl 



e f2 



F -l T (k 1+ k 2 ) 



N 

n 

n=l 



fe/3 



r 3 



u2 



i'ii 3 



In general, for the mth SU the conditional PMF is 
p(k m |ki,k 2 ,...,k m _i) 



F f - v; : ' fc 7j 



F 



1 T (E^ 

PS 



n 



Ei=i fc 

knm. 



"J 



Using the chain rule and factorization of PMFs, the joint 
PMF for SUs is expressed as 

m 

p(ki,k 2 ,...,k m ) = JJp(k r |k r _i,k r -2).--)ki) p(k x ). 

r=2 

Finally, the marginal PMF of the mth SU with multiple 
N PUs can be obtained. Based on the evaluations above, 
it is straightforward to obtain the expected value of k nm . 
Therefore, it is omitted for brevity. 
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